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(54) [Title of the Invention] Computer-assisted 
Simulation Method Obtaining an Electromagnetic Field 
of an Object 

(57) [Abstract] 

5 [Object] The electromagnetic field of an object is 
obtained with a computer, so that" the electromagnetic 
compatibility of the object is ensured before the object 
is implemented, 

[Means for Solving the Problem] Global multipole 

10 expansion which represents the operations of an electric 
charge and an electric current to a distant point in 
each partial region in the form of multipole expansion, 
and local multipole expansion which represents the 
operations of an electric charge and an electric current 

15 to a point in each partial region are performed, and 
the global multipole expansion and the local multipole 
expansion are superposed, whereby the magnetic field 
of an object is obtained, 
[What is claimed is] 

20 [claim 1] A computer-assisted simulation method 
obtaining an electromagnetic field of an object including 
a plurality of regions, a plurality of charges and a 
plurality of currents, comprising the steps of: 
a) performing global multipole expansion which 

25 represents operations of the plurality of charges and 
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the plurality of currents to a distant point in each 
of the plurality of partial regions in a f orm of multipole 
expansion, and local multipole expansion which 
represents operations of the plurality of charges and 
5 the plurality of currents at a point in each of the 
plurality of partial regions in a form of multipole 
expansion; and 

b) obtaining the electromagnetic field of the object 
by superposing the global multipole expansion and the 
10 local multipole expansion in the plurality of partial 
regions . 

[claim 2] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 

15 to claim 1, wherein 

the following steps obtaining an electromagnetic 
field of an object are repeatedly executed up to a given 
small error scale, 1=0 satisfied as an initial condition, 
the steps being: 

20 a) calculating the global multipole expansion by 
using a global multipole coefficient based on C^=GI, 
being a vector having a global multipole coefficient 
in a plurality of partial regions, I being a vector 
representing a given current distribution, and G being 

25 a matrix determining a global multipole coefficient in 
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one corresponding partial region among the plurality 
of partial regions in the given current distribution 
I; 

b) calculating the local multipole expansion by using 
a local multipole coefficient based on c^=Tc^, c} being 
a vector having the local multipole coefficient in the 
plurality of partial regions, and T being a moving matrix 
used to form local multipole by combining global 
multipoles; and 

c) obtaining an electromagnetic field f rom ZI=Z ' +Lc^, 
Z being an impedance matrix, Z' beingpart of the impedance 
matrix Z representing a connection between the plurality 
of partial regions, and L being a matrix evaluating a 
local multipole coefficient. 

[claim 3] The computer-assisted simulation method 
obtaining an electromagnetic field of an object according 
to claim 1 or 2, wherein 

the plurality of partial regions are identical in 

size • 

[claim 4] The computer-assisted simulation method 
obtaining an electromagnetic filed of an ob j ect according 
to claim 1 or 2, wherein 

a size of a partial region is formed in proportion 
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to a distance from an observer region* 

[claim 5] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 
to any one of claims 1 to 4, wherein 

one corresponding partial region among the 
plurality of partial regions is respectively allocated 
to one zone having a uniform physical characteristic, 

[claim 6] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 
to any one of claims 1 to 5, wherein 

one corresponding partial region among the 
plurality of partial regions is partitioned into 8 zones 
at a maximum. 

[claim 7] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 
to any one of claims 2 to 6, wherein 

an element which is a constituent portion of a 
partial region of an object and has an impedance is 
directly considered as an impedance directly in a matrix 
Z' . 

[claim 8] The computer-assisted simulation method 
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obtaining an electromagnetic field of an object according 
to any one of claims 1 to 7, wherein 

an electromagnetic filed is obtained for a given 
frequency. 

5 

[claim 9] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 
to claim 8, wherein: 

a given frequency is determined from a minimum or 

10 a maximum frequency; and 

the computer-assisted simulationmethodobtaining 
an electromagnetic field of an object is started from 
the minimum frequency and is continued to the maximum 
frequency in a given step width, and an electronic field 

15 is obtained at each time. 

[claim 10] The computer-assisted simulation method 
obtaining an electromagnetic field of anobj ect according 
to claim 8, wherein: 

20 a given frequency is determined from a minimum or 

a maximum frequency; and 

the computer-assisted simulationmethodobtaining 
an electromagnetic field of an object is started from 
the minimum frequency and is continued to the maximum 

25 frequency in a given step width, and an electronic field 
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is obtained at each time. 

[claim 11] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof anobject according 
5 to claim 8, wherein: 

a given frequency is determined from a minimum or 
a maximum frequency; and 

the computer-assisted simulationmethodobtaining 
an electromagnetic field of an object is started from 
10 a frequency between the minimum frequency and the maximum 
frequency and is continued to the maximum frequency or 
the minimum frequency in a given step width, and an 
electronic field is obtained at each time. 

15 [claim 12] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an object according 
to any one of claims 1 through 11, wherein 

the global multipole expansion is performed not 
in accordance with a basis function but for each element, 

20 so that stability at a low frequency is ensured. 

[claim 13] The computer-assisted simulation method 
obtaining an electromagnetic f ieldof an obj ect according 
to any one of claims 1 through 12, wherein 
25 electromagnetic compatibility (EMC) of an object 
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(54) COMPUTER-AIDED SIMULATION METHOD FOR 
OBTAINING ELECTROMAGNETIC FIELD OF 
OBJECT 

(57) Abstract: 

PROBLEM TO BE SOLVED: To guarantee 
electromagnetic environmental compatibility by obtaining 
the electromagnetic field of an object with the 
superposition of the global multipolar deployment 
Indicating the action of electric charges and currents 
to a long-distance point via multipolar deployment and 
the local multipolar deployment Indicating the action of 
electric charges and currents to a point via multipolar 
deployment. 

SOLUTION: The electromagnetic field of an object 
containing multiple partial regions and containing 
multiple electric charges and currents is obtained by 
the computer-aided simulation method. Global multipolar 
deployment and local multipolar deployment are 
implemented in the partial regions, the global 
multipolar deployment indicates the action of electric 
charges and currents to a long-distance point via 
multipolar deployment, and likewise the local multipolar 
deployment indicates the action of electric charges and 
currents to a point in multiple partial regions in the 
form of multipolar deployment. The global multipolar 
deployment Is implemented for the investigation object 
sub-divided into partial regions in step la. The local 



multipolar deployment is implemented for the object 
sub-divided into partial regions in step lb, and the 
electromagnetic field of the object is finally obtained 
via the superposition of the multipolar deployments in 
step 1c. 
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[0 0 6 9] f) -eti^tKOffiSBBI^^r^EEiK^ttS 
[0 0 7 0] 

i^*^#i!)-?5^lCfcV ^•C liJ^^4^E tan= 0 tt?:^?^ ( 3 
0) j: y^$tlS«S^<Z)?^T'^Stl-5o 

[0 0 7 1] 

^1 0] 

(30) 



[0 0 7 4] reSlC^bT 
(i coA+V*) tan=Eitan 
fefe'U Alir^(3 2a) lCj:y^3tlS3J<7^>^'-V 

;k <^m9t^ (3 2b) icj:y^snsrf<7^>2/^;i/T' 

[0 0 7 5] 

) mi 1] 



s 

-ihR 



(32b) 



^(31) tmm^i^yri.m\z.umm^Kt>fz^n^m^f^^t^ ( 

3 3) 

+ V<p).dr « J • dr mitTn=i.2 M (33) 

Cm 

f^^iLm«=l. 2 M-C*4o 



[0 0 7 6] -C0^lCV*lCgaLT«to4xS^S:XSl ^>hlC*fLT?^^ (3 4) im^tl^o 

i5<DmMM^(Dif^'ry'y^)\^\^^'')^^^^^ to 0 7 7] 



(9) 



11- 38059 



15 



16 



(34) 



/;/iLm=l, 2 M-C*4o 

[0 0 7 8] M-^<Dmm&m^i^ hJld^alJSC^^^^n&l^ 

[0 0 7 9] 1 r3(7)^rj^li, ^(2 5) 

(Dm^^m^^^WL. rolled: Ut1^*^:^ifSlCfSm 10 

[0 0 8 0] fm^m^o)^^ * 



(3 2 a) JSt0 (3 2 b) ©J^T^^S/^^WCfeVNT^ 
(2 9 a) 2fctJ^^ (2 9 b) «:^bT?^^ (3 5 a) 
JSOJ^ (3 5 b) iCjty^StlSBS^^^^ejtlSo 
[00 8 1] 
mi 3] 



N -ikR 



n=l B„ 



<p(r) 



fc/SLSn : =s„iUS,;2'^*>ao 

[0 0 8 2] m.-oX. *feimin=an^^^'^S^- 
(34) tt?^^ (3 6) ©?^(DVhU^':;^:^S^i:L/ 

[0083] Z • I=U (36) 

fc^cbvhy^'xz= [Znin] ii-f>i^-^f>;^-7h'; 

^y^.T'feS, ■|pri!iS[35:e3liZ»®!^Si|iSl= [In] «:'fiS 
U= [UJ lC||^bTlr^<5*^e>T'S)So IS^lCiS^Vh 



2 Vg • ifrn(l") «la' 



(35a) 



135b) 



5?J 



85 [0 0 8 4] #is?)S^-^>hIn»:&@^ O 6) 

[0 0 8 5] ffiSn%2oCD®xl/:jl>KSnlWSn2 
tC^^SIb^i^tC^ (3 4) 3!)^e)VhU^'X^ll/^>h 
(VhUi^Xl^) ZmnlC?t*UT?i*C^ (3 7) *^#e>*l 

[00 8 6] 
BS14] 



= ia. (ftii + aji) ^ - (Al2 + ^22) + 

fc/iL^f^^-^-JVApqlcWLT;**:* (3 8 a) *«SE&*>> sJf-rVv-i-JV^ 
lzW-UT»i*« (3 8 b) i«tfa!:-3o 
ikKq 



pq' 



<Ppq 



Sap 



Bap 



<38a) 



(3eb> 



[0 0 8 7] u^^y\f~'S(y7.t.^^^ti^'rhV^ 

>^x > h Z nn»llm«)X 1/ > hM^CDH n (DX 1/ 
^>hM«)i^=&^t„ Ei4^CfeV^Ty-^Xl/^>h 

€>o 'WIA«A2il*, Sn2t-<ky Snil«)*,#.lC:^$nfe 

[0 0 8 8] >f >l^-^>;^>VhU^:^ [Zmn] ^tiM- 



50 



7) ICiyM^^^ZnmSrS'JM^^o^lCfl^fffi-ra^i: 
{CfeSo U*^LH^?^xL/;jl>h®#^lc;lilO(Z)x^ 

il#^e>tl-g)o ZlCDfciblC^ (3 8 a) JtCJ^^ (3 8 
b ) IC J: SgP^^^jb^^tL^n 2 ocDH^J^x 1/ > h 

U 7>x > h tcinl^^ tlSo 



(10) 



11- 38059 



17 



18 



[0 0 8 91 (3 6) (DWit^ (3 4) tCJ: 

Ur^iJit (3 9) (Z)J:e)tCftit3tlSo * 



[009 0] 
BSl 5] 



C39) 



[0 0 9 1] J^Ki; h)\yt.LXTf.i:^l^UWS[»ltL 

xEomi■&^M^'mmzn\>^xE^mk^ (40) \z 

[0 0 9 2] Ei (r) =Eo • e-ikr (4 0) 
[0093] Etan=zs' K 1) 

wzimt^^ \^i>^\M^mmz.i^^^xmm^-&^ 

Cm 



10 



Cm 



[0 0 9 7] fcfcXiC; (4 3) cD;63a®:^^^S:MlES 
[0 0 9 8] 125 a^C^3V^■C:^>^S^lfe^>tr-^S^>^ 



JzsK 
Cm 



dsr = Z In j" ZaVn 
a Cm 



[0 10 0] ^^eoTlo®H^?^xl/J^>h^uilS« 

(3 6) CDif(Z)9o*-t'CD"v'hU^Xfi8;4J-*'^ffiiES 
lo(Z)^Xl/^>hSUjiSlil 6*T*®vh'J 

[0 10 1] J:yfa#icii, :iyif^%yVttAlo'^m^ 



>>^zs&<gv>jg?sgg5«cfcv\T3a<emtcr;j(^ (42) tc 

[0094] z s^rsp/ld (4 2) 
[0 0 9 5] iTLV>J^^^ (^(4 1) $:#M) ICcfc 

(3 0) ®«*^#i'tlfc9-=fo?^^ (4 3) ic^'fb-r 

[0 0 9 6] 
[IKl 6] 

(43) 

MiES (-r^fct)^ (4 3) (O* 
j2) lC*tbT?^^ (4 4) 
[0 0 9 9] 
OKI 7] 

(44) 

m^-c^i:^ig^ (3 6) (^ifcDfefc' i 

Ay\f.-9yy^-L\^M-t^^^ 

(I^:iLvV\ZM\uXX>f^ (4 5) 
[0 10 2] 

m\ 8] 



30 



11=1 

[0 10 3] f=.tlV^ (3 7) C!)*CDi®Zmni:^ (3 40O-53i:*^liT'feSc 

9) <z)t^I(z>fitUIni:^*^«)e>^l•rv^■5= Mi&n^Sr^saic [o i 0 4] 

^lir-rsi:, z:(3!)#*II«5:, (4 6) {Cj:U*3n OKI 9] 

[0 10 5] ^if-fbstifc^Eiii 



(45) 



(46) 

(3 6) o^tctii-omrimm^^^-^^^ 

[0 10 6] 
[gS[20] 



19 



(11) 
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N 



(47) 



n=l 



- Tjl <'*8) 



[0 10 7] mf^<^^^ 

mmi)'^mmx'&?>m^i^i-i-^ i3W^ (3 6) ^m< 

[0 10 8] -(D*^Stifc^lc<};y^^^*<D^4)- 
[0 10 9] EMCe/^3.l/-e/3>a):J>S^T*r®J;e) 

[0 110] m^ommmm 



da' 



-ikR 



da' 



s 

[0 113] Z:tie>O^l4R=0X«r=r' ^cev^ 
5= 



io*>r>i^-^>^-7hu^:^&fm-r'5fc«3tc^ (3 8 
a) (3 8 b) <m^^mm\z.mmt^z.iiii'^'Mm: 

[0 1 1 1] ^ (3 8 a) 25.t>* (3 8 b) (Z)t|JlC?^^ 
(4 9 a) S.?>* (4 9 b) ICJ: y^$tlS^>f ::^CD^Ii)' 
20 fm&t&o 
[0 112] 
B&21] 

(49a) 
(49b) 

[0 1 1 4] ^ (4 9 a) a)a^^•^CfcV^TV^<o;6^0^ 
30 ?^«)^icr^ (5 0 a) *^e>ti.s= 
[0 115] 
QK22] 



21 



(12) 



e^m^ 1 1 - 3 8 0 5 9 

22 



-ikR rT^I^ - 1 
S J ' 

-ikR _ 



I02 



IG3 



IG4 



Ai^tr-Ji. ^ (SOa) Ci3V>Tti » "C^^nx^.^io 
^ (4 9 b) £?3|S^t;*JV^-ClWCge?^$n/;^ (5 0 b) iim^tl^o 



IJ ^a- oJJ da. . JJ i da- <SOb, 



S ^ ^ 



1G2 



[0116] mm^ I G3S;t>' I G^^mmmizmstt * [on?] r-^ 0 iz^its^ipsfiti: lx i G2®^t« 

liX l/^ > h^®^^iC<D*-eitiggfi^T'*> y ffJRSS tlT V ^ [0118] 



lim — 



-IkR 



-tk 



-ikR _ . 



lim 
R-^0 



(51) 



(52) 



[0 119] m9mm 

m^ti^:^m^ (3 6) i:wwmi^^^m<i5mD 

^ (3 6) 7!)^e,(Dft?<^^»*@^^±-fl$S5^- 
(5 3) icj;y^§nSo 

[0 120] Ax = b (53) 
•to iaie(«¥x(n')cD^$tlfc?iJT' 

[0 12 1] mm^-h"^^ ismz.^^ht\.f~tm:=i7.v 



40 [0 12 2] ijW^^m<r=.^\zm^Mt. m^<D 

Emmz^^LX 1 U l oRi:(Z)v 

v^vhU4'>^.®^^c^i3a)fe«)^ciiS^io (n2) # 
z 50 S|;fc^««iJe^T, mM<^S^:^7^^v>^(Z)iSS:T'gSfc'W 



(13) 
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[0 12 3] CGS-r^t)*>^t^S!M^54iScDi^^;^{cli, 

Youcef SaasS.martin H. Schultz*" GMRES: A Gener 
alized Minimal Residual Algorithm for Solving Nons 
ymmetric Linear Systems" (SIAM J. Sci. Stat. Comp. 

e)lC-rSfc«)tC, ^ (5 3) f)^^imhXrf^ (54) 
tC ^ y tl<5S!ls^ r (m) /)^H^t!Cx (m) a)7ti?)<^^ 



24 



(5 5) lCj:Ua$tlSl^(Z)^e(ni)«, iE?f| 

[0 12 6] e (ii') = x-x(i") (5 5) 
irj^ r (m) i^m(DWMX'^^. "T^t^-b^-^l > 

■r^WBtV&^o ^ (5 5) (5 4) (©cfilCf^ 
U ^ (5 3) $:#Ji;-rSi:?^^ (5 6) ICjcU^Stl 

[0 12 7] r('») = Ae(n) (5 6) 

J II r (>») 11 m^oMi t m7L\t 1 

m:^^ (5 7) ICJcy^StlS^P^A^MSnSh 
[0 12 8] 



20 



[0 12 5] r(iii) = b-Ax('n) (5 4) * 

II r (n) II / II r (0) II = II r ("i) II / II b IK £ (5 7) 

[0 12 9] ^T^?3^J^^cPM■rs*T'^CiK«^M«:^ 

2, . . . , MT'fe^c 
[0 13 0] ^Un^'-r-f i/3— 



[0 13 1] ^m(D-:f^)ny7'^i/^—y^h&MC0-f 
(5 3) ij^fi^. t^m-^titcU^^JhizMLXik^ 

(5 8) *^etlS= 
[0 13 2] 



b = (SB) 



[0 13 3] '7\^V^7M-Ut'fV::iy'f^ty3-^t.^ 
1 T*fe 'A t^oTiE5l^ft¥lifefe' 1 lHl<Z)H^a)^lcSElc 

ti, V h U ^ A ^1 V h U ^ :z.M}:Lim\&X':^ti\,t'fV 
[0 1 3 4] ^ (5 8) S:^M■rs^^c|ib"^iS^ 

A"ttSii^®aft*^e>5tfiJlciiW5tS^tcii|ti[3ti^ 

a U"^ X ts to & 4 

[0 13 7] ^ (5 9) ?:^-r-5#^. :^«bS^x$: 

;i/«x =M-1 X * Srfflv NTl-f^UJfcl-ttiliJfc e);&v ^„ 
[0138] l*A " (DmmU 50 



[0 13 5] &m(D^V:ny7'^i/B-y'^tmmiz:^ 
40 {)l!l©:/U=l>-r'f i>3->W^ (5 9) ®cfce){C^f 

[0 13 6] 
[li$C25] 



(59) 



M-^ay^iti\m(Di^mzmbxm^xsb^o b*>u, 
4 ) aM\^i:w^<DRmt LT-eo* ^arnxmrn-t^ 



25 

=Lyf(Dm^K\i^ (5 4) <Dm'0i\t-m»i^(Dm^\Z.\t 

\Z\im\\Z\t^ ^ (5 4) ICfeVNT^^^tlfcftJ^C^* 

&(<») = b-^«(»> ° n-M*) 

[0 14 0] z:<D#^iciiT h U i'^M-itcj: 

[0 14 1] •Yntf:/Un>f^-<i/a— >^ 
1 o(D oIteJ&:/U n >f^>f J/ a =.> ^^li^:?:^ (61) fc 

[0 14 2] i = j©#^lCl*Mij = Aij 
-e®flfi<Z)«^lCliMi j =0 (6 1) 

[0 14 3] lo^M^S^tftJ^^U^^T^-f i^3-> 

F= {1, 2 N} 

[0 144] 

w 

[0 14 5] ^CJ:U*^#^C2V^^C^T^^'5gI5^-^Fw 

ricv ^■C'■=r V^)^ 7.ymif(^ (62) ic J: y e>4if>. 

[0 14 6] i. jeFw(7)#^lCl*Mij = Aij 
-ecJDftfiCJD^lCliMij =0 (62) 

[0147] 1 o(Z)i^lc^^?^iR^®#fflli> yrv 

U -^/5i6t<Z)?^{Cy^.7^^v •;^#^c-^^:/^t|■t-S 3 tc J: 
IC:^S tlfeit)?«s«^ftJl-li^ 1 iac^^ ® 1 / 2 & 

^cfe^.^Tlo;e^'t^i*D-re= 
[0 14 8] ^Jicii, ^(D;5m:7'n-fe;^.©|trtcgfi 



I jj^ipii- 3 8059 
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(6 0) 

[0 13 9] 
[SS2 6] 



(60) 

tj d7>t*^-7 h U XM-l^^figStiSo 
[0 14 9] EI6lCfeV\Tmi£§i&*®*<*^«^:J>S'J*^ 

iM^'fb>^7^>y::^2. 3S.t>*4^c^3v^•C5^$tlTv^5„ 
[0 15 0] •Ynif:/U:a>7^^' i/3— >^^«)7iffll*> 

SSlSKiaic J: y T'#Sfe't-tiE5Stc:g^$tiS <fc o IC 
lciigp^IJ^li-eti-?ti®MJi!cl*lfct)l4^3 OQ^M 
[0 15 1] FFS^ (ffi5iM?adg5c^7^'yi^>iO 

[0 15 2] ^^6^^&^1--5#^lCtt^3XhS: 

[0 15 3] 1 oC!)3g^S^lJ:^^T'«X^- hfli 

IJ, «^(Z)3o(D«?^i' N;bxi-3, xi_2;5tt>*xi-i*^ 
e>>^^t-hfflxi(0)&2^^fSr^ (6 3) *^#^E 
50 f -So 



(15) 



27 



[0154] * 
xi(0) = xi-3- 3x1-2+3x1-1 (6 3) 

[0 15 5] Z(Diy:iL-^tti!i&(DW^r^O-^il(D 

#^C!)fe«)a)T;i/=fUXA*V:?C^ (Alg-1) ICiU 
[0 15 6] f : =fmin 

M-1 : =A (xiiAcoais^) ^tmt^ 

X : =M-1 b 

f : =f +Af 

>(.^l-hMx'(0)S:ffl^t-S 

AM-1 X ' = b ' ^BL^m\zm< 

X : =M~lx* 

mi (Alg-1) 

-Ha®! mii^J A««*J<SV ^JSM f mint'X - h S 

U#?>tl<l>= AC[)|-t«3nfe»®?^^®^«?^V^T% 

^M-iilbTjgfflStl-§.= ZLCOfciblCMMf li^^^v 30 
^^tCfn,in*^e>fniax*-ei««)e5tU ^fefflB®^ 

CDlfetCM^i^^x : =M-lx'$:/rbT^«)e5tv53il*^ 
[0 15 7] m^(Dtm:=i^ hltFF S&lCfcVNT^g 



i^F^^ 11- 38059 
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[0 15 8] m^<^M\m.ii'^F¥sm\ziQ\,^xvmtiti 

}b^y--fTy7.\mf\^^r)V^^)XK (Alg-1) \z 
^StlT V^■5)<fc \Z f min^'^j f max^® 

[0 15 9] 3iic-:;^i;3>-r-<S/3-t'M-l&¥i^MifigS: 

\mmiKi^ix,tS^^ f max^ f max*^^ f min'^Xli* 
l^felC f ctl^h f Biin^?^*'^"^' f c*^e> f max^ilMS 

-5 r. i: j&^tgT'feS z. {cfeSo 
[0 16 0] ji^;KT->i/-v;i/CO#^iilM 

^^w^oMWiit. ffmoymk^^(Dif^7-yt/^)\^'BL 

[0 16 1] :$:5l^lcevNT^$tiS^<5*^#lt^fe 

[0 16 2] z.t\.<Dfmh\^x. mm(D^-^E(D^mM 

[0 16 3] 3mm(D^k^ (6 4) lCj:U^$tl€>|Bir^ 
•5, 

[0 164] (V2+k2) •u = 0 (64) 

[0 16 5] :itih(m^m\'^x^<DmLui:imm. 

[0 16 6] ^;i/Ad-N;i/^y^^^ (6 4) ^ms<D^m 

©r, 0, a&MV^T^t^lCliii^ (6 5) 
[0 16 7] 

Ba28] 



(65) 



[0 16 8] ^^^■;?5g^&fi?<^-*?3l-?^^ (6 6) IZ 
[0 169] 

u = Ur (r) -Ue (0) *Ua (a) (6 6) 



50 



m(Dmzimti>o ^(6 6) S:^ (6 5) ©tftcf^ 
r2 s i n^d /u 

iim%-t^:r.tizj:*)i!^ (6 7) *^#ejtlSo 
[0 17 0] 



m29] 



29 



(16) 
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Bin's d r^3fc^ 4. Sill ^fBin»^l + Bin* « = 0 

(67) 

Uo da' 

::tLUJ:»)* (6 7) iPh. o*«'b»fi:LTv>«ratM: <6 9) *<»f,Kao 
d7 ^ dr^ Oft Bin » d» ^ sin* 8 



0 

(69) 



[0 17 1] e^(mm^^m-^t=-i^i^^ ^ (6 * [o 1 7 2] 
9 ) ma^mmo^mif^mm uc j: uriJc^ ( 7 0 ) bs3 0 1 



Oft Bin » d» V dfi > Bin^ » 

^ (7 1 ) ij^tsna r ic^fi^-t4tKa-:^^*^*= 

J_ A - 1(1 + 1) + k V = 0 

Or te dr 

[0 17 3] ^6^^iZ^*)^ (6 6) *^e)^;^A5J^ 58 [0 17 4] 

)iy:^m^ii. s-DCD-m^ism^ (72a), (72 msi] 

b) S.t>* (7 2 c) ICl^StlSo 

A(^^).[(.r)=-Ki.i)]u. = o 



(71> 



(72a> 



Bins «A 



da' 



[0 11 5] :^m^ a 2 c) it. ^<DmC^mti> t^[0 17 8] 

SM^>;^SiCT'fey, (73) tCj:U^3tiSj»5: Ue.im-Pl"' (cose) .Qi"! (cose) (74) 

[0 17 6] Ua.m ~ eim" (7 3) 40 VNTiim^.$:^U tJeoT^^lJlCli^^tl^o 1^<D 

2im^Ua (a + 2^) =Ua (a) J^'J^ ml*^ miiBiRPl'" (cMse) T'feU, fefe'U 1. mGNo 

[0 1 7 7] (7 2 b) ii;i/i/>>K;i/:&@^®io [0179] 

[0 18 0] z.mz^^jmtm'^n^, pi® (cose) =pi (cose) (7 6) 

[0181] 50 R-Cfm> 1 



(17) l^m^ll- 3 8059 

31 ^2 

Pim (COS0) =0 (7 7) * n*^^«)/iz;i/0§IK$:*-r€,^'y-fe;i/1»4^:&S^T«> 

mmz^^m^ti^o [0 18 3] 

[0 18 2] PqO (cose) =1 Ur.i~Ji (kr) . yi (kr) (7 8) 

PiO Cms 9) =cose z:tie>0D'<^yfe;i/5iKl^ ®E^(Z)Rtt^'v-fe;bBiiSJi 

p i (cose) =-sine (X) . Yl (X) K^mm^ (7 9) S:/M.TK54 

P2O (cose) =3/2cos2e- 1/2 #lte,tU 15J^ISIC j; 
P2I (cose) =-3sinecose [0 184] 

P22 (cose) =3sin2e ^3 3] 

:^m. (7 2 a) ljU »?i:bTr^^ (7 8) icj:y^$*io 

jl(x) " • Jl+0.5(X) , yiW = • *L+05(X) (79) 

[0 18 5] r . ^ ^ ^ 

0186] jo=sin (kr) /k r 3^K^y*;Hiaf!t$:ffl*^S^i:lC (80a) 

(kr) /kr S^y^* (8 0b) lCj:y^S*l5lllS&tKII2S(^^N 

j i=sin (k r) / (k r) 2-cos (k r) /k r >^;bBBI5C*^?5n€.o 

yi (k r) cos (k r) / (k r) 2-sin (k r) 5« [0 18 7] 

hl(l)(kr) = ji(kr)+i • yi(kr) (80a) 
hi(2)(kr) = jl(kr)-i • yi(kr) (80b) 
±iE®B3|jtli^^u®^JilCKLT, ggiSJl (k 20T!!rl»hi(l) (k r). hi(2)(k r)«r^0{C*fbTi|#MlC 
r) . yi (kr) (C^T'feSo BSii5Cu(Z)^lCft# «c«.*^e>T*fe^. 

UT^^^y •fe;i^gltXli«^>>^;«/B8l»&^^^*'2> [0188] ^©cD|lI^^:^;I/A;i^;^^y:&S^®#«)5 

i±B9t!cj 1 (k r) (D^^mm-t^rd»r'^X'&^. [O l 8 9] 
mtjt^\mmyi (kr) jmx>^^'citm^^y^^^^ ^ 

uiai=bl(kr) • Pi I"! (cose) • eima (81) 
fefe'bbi(kr)timir^<i[)^'<>yfe;i/BiSaXl*«>'^> (8 2) {Cj:y^^*lS»SM*^e>tlS„ 

>!r;bB9l!C(7)*^iK^S:^U lGNo^tt>*meZT'fe [0 19 1] 

S„ [|S3 4] 

[0 1 9 0] ^^uiC^LT^ (7 7) ^^mLriX^30 



eo 1 
1 = 0 in"-l 



1 = 0 m"-l 



<B2) 



[0 19 2] tc:t£LmiCin>it^n:^$>^. OmS 1 C0#^IC1±. «®BBm5:?k^ (8 3) ICJ: U^f 

[0193] mWm&L ^ i: A^T'SSo 

it:; (8 2) CD^SlC^Sc#-rS«Pl I " I (cose) &t;e [0 19 4] 

ina^rf^fli^^L-^Sil. V>t>»i>S«Sg9IS5CYl'»(9, a)*^ 40 [183 5] 



gina (83) 



[0 19 5] laKOOmmmLmk^ (8 4) (D^vffi 
{CiU^^tlS, [0 19 7] 

[0 1 9 6] Yrin(e, a) = (-l)«'Yim(e, a)* [3K3 6] 

^nic±yr>cic (8 5) tcj:y^$n-5jE«i£^^w* 

{», a) • T^'(8, a)* «to = «11' • «imi' 



n 



(18) 
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[0 19 8] fcfc'u Q},tMim(D0m. sit^m (8 

[0 19 9] 1 = 1' (Dm^iZlidir = 1 

^<m(Dm^izit8ir=o (8 6) 

3^®g8iaii ^"^'^ yr^if ify^l JScC^mlzm UT^jE 



*tt®io*^e)ti<g), t-'fct)*.ffit<o^^llSISg(e, 

a)$:MKISYlH(e, a)S:ffiVNT??:^ (8 7) (D^ 
[0 2 0 0] 

ms7] 



1=0 m=-l 

m:CiJi^ (8 5) (88) «cx>>*fte>n«„ 



(87) 



(80) 



ft 



[0 2 0 1] mor^*^ 



2r^C*T*<3[)OTBi»li«^l- 85 [0 2 0 2] 

J« 0K3 8] 



1 



^1 c lii . Bin d ' cos S • e 



-ia 



^ V8JI 



yO = J— • COB « 



= -,|— • Bin » • COB » • e 



la 



Bin' 



(8 9) {cj:y^snsyu->si5c mz^<'mm-r^m^ (90) *^#p,tiSo 

G(r, r) = e-ikR/R (89) [0204] 
(f-fe' b R = II r - r • II T'fe<g>) -S(^)^* 3 9 ] 



,-lJtR 



1=0 in=-l 

(90) 



[0 2 0 5] ^wMm 

^ (9 0) &fflv^T3ij^3i<7=->i/^;bA:at»(o#'^ 
3SS:#m-r-5z:i:/)^BriiT'fes= zLCDtc^mz^ (90) 

(14 a) ^lt>' (14 b) (Z)cf l-f^AU ^t.M 

'(pi«fcl:e>«r = r' t^^-^smtstm^imiz^ (9 

0) *^e>l^3tlTV^S*^eJT'$>-5= h;i/5K^>i/ 

Y;i/A<3DilBlli, ^ (9 0) ^rSocDiCi^SJ^m^ 

Ax, Ayjscx:^Az.iz%mizmmi-^^tiz^*mt>n 



[0 2 0 6] |^#<^?jtt)^y-:x.;^,®aetC'«K#LTM 

V^) mBm^<^t:^^0>^t. (^DJ:^J:*CgV^) 3^ 
50 ®||li(Z)fefcC0i^i;S:^1-. 
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[0 2 0 7] m&ii<^yt/^J\^0)^t 

(9 1 a) lCj;y^S4x-2>„ 



(19) 

* [0 2 0 8] 
* [85(4 0] 



1 1 - 3 8 0 5 9 
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loO tn=-l 



(9la) 



*0 1-0 m=-l 

° (91b) 
Q 

J^^^ i^^r)VAjm * lC*f bT«k^ ( 9 2 a ) 

[0 2 10] E19<Db^C(5&5t$tlTV^S^> -tJaitot) 
S^#.^,5b^v-:5^,^.t-lt bTMm^,(^i£< ICfig-tS 20 [0211] 



/^jeL±ES:<r»»>ff30-*/v#e^l^raLT*^ (9 2 b) *«ffifiK'o 
,1 = J/J J(r') 4^^^:-) o ) dV 



(92a) 



(92b) 



[0 2 12] >^oy^i^^65^^CO|gfflT'«, 
JSgg (9 1 a) ^tt>* (9 2 a) itmmoMKDTMfDWt'tJ 

n^mhm<D^mmmt>mz (l+d 2(^^s:#t? 



5. tm.it. d/D*vJ^s^m^i/^sv^^^cJ:^Ja^^^c 



5^>2/^;i/A, *(D#^M®fe*i)lC^|JT'4 (L + 

1) 2ayy^iJ^m^tM!mt•^i^?>o 

[0 2 13] 1 (Z)ffitcmbTfi:t)ail«e^^^i±iS 
[0 2 14] 1 = 0 : 

1 = 1 

1 = 2 
1 = 3 
1 = 4 



2 
4 

8_ 
1 6 



- > '2>fc«?)lC^JlCtt L ^ k d 

[0 2 15] ^jmMm<D^mmm 

[0 2 16] •*}-^^mi}^^n^ti^'^ti'en<7>mmi^ 
50 mmtmion^. ■t5fc*)t.g*SFw®#^cD>f>-r^y 



(20) am^ll- 3 80 59 
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JC0gbT^ft)tlS„ [0 2 18] 

[0 2 17] »£M!KW-n (2 9 a) ^t?>' (2 9 BS4 2] 

b) $:#m) lCj;S®jS^^®^Jg«:^ (9 la. b) * 



W = -u*o i Z 4n.,w (8. a)' 

1=0 mo-1 



(93a) 



1 = 0 m=-l 

f.fil_jjg^n!p<r>^n-/<A,^mS>mi»A (9 3 b) 

(93b) 



[0219] H{C«tJzaG)f^lC J: y ffl^ ©SPf^lclCfe }« [ 0 2 2 0 ] 

i:rsa-*;i/#=figiBg*v^ (9 4 a) lCj:t;^Sti ^4 3] 

5« 



■5 



= -ik^o i i: ii()cp) vr(s, «r 

'^0 1-0 Bl'-l 

/:«LJiE3W)'«'<iCin-*JV#g*^«4^ (9 4 b) 
ai^« = 2 In iJ Vn(r') 4'^(l^) vr(S'. «') da' 

<94b> 

nEFw\II„ Sn 



[022 1] ^©r, eR-Cfoc(D±(D9-^ ^ 

ziiife^fo (9 4 b) 0EtJCDn-*;i/#^^ 40 

[0 2 2 2] EMCi/^a.U-i/a^COflflicWcD^ef 

{^wmmL. mm) \tm^\\zm^-&\<^mm. 

m,Zti^^<DX\ ^ (9 3 a) -BtXf- (9 4 a) 

yn->'t;i/#OTM (9 3 a) A^g^I^o^gPT'fNi 

s*i-5#^ic«, 3J!?cf^r4®5i:^{*2K©gi54^iJac{cfe 

50 



[0 2 2 3] 
B55C44] 

d ^ X 

[0224] ^\,^mm^ziQ\'i^^'^'mmm 

[0 2 2 5] Ell o\tz.(D^^m^tmw)\^—:f^m 

icLT^-fo zico?|i(!}c;i/--;''a):^:^93K5^>J^-^;i/tt4 
0(55^^0 -/'<;i/#1iJllB0ic J: u^stis (EiiocD 
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(2 9 b) #fis) (Dmn 

[0 2 2 6] :i<D®il«IIBK<^)#^iSIW*'^f^tiS*frf- 
micsv >tc^^ >-fe;i/-r-6a^S:-^;^S zi ic J: U 35c 

[0 2 2 7] (FMM) 

okhlinOH^i" RapidSolution of Integral Equations 
of Classical Potential Theory. " (Journal of Comp 
utational Physics, ^6 0#, HI 8 iM-^mz 0 7 
1 9 8 5^^) ^f-f-^-5^V^*c#^ 

<Dmm^ji.^^izi5^^X" ^yf)\^" S*u ?:^:v^T'■9->:/ so 

m<m±itmmzx'^f£^^(Dxmmt.^j^^:r.- u r 

%228] rz:-C'*^-tS*^l*4^l- 

[0 2 2 9] t.wmi3^o^\ts mi^^tiEMcm 

[0 2 3 0] S(fijEcDH^ft??*$:> S©^<»^*^J^ 
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•^h^)i!;7s' V)ym>\^x\-i.W$i7.'T^v:ftt!^o 

(N2) #i6/JMSfo?M^«fc 'J ^f-^^ V ^?t^T'^f -r S r il 
[0 2 3 1] m^#.«?^^ (9 5) IC J; U^$tlS?^T' 
[0 2 3 2] Z = Z' +LTG (9 5) 

^-f£u z' ^d:vhu^;7>^Z(z>^v^^c||^@^■ts^l5^I^i 

®r^<Dii@^S:^-r§l5^5:a^U LTGliv h U i^XZ 

;fct)t)iGFw. jeF\Nw§:5^-to 
[0 2 3 3] |tfjZB(Z>^S^[:^* ({SVNMMt-fcltS^^ 

XZ' tt?^^ (9 6) lCJ:y*«?)btl-g>. 
[0 2 3 4] 

iGFw, jeNw®#^tciiz' ij=Zij 

^CDffi^^lCttZ' ij=0 (9 6) 

So 

[0 2 3 5] 01 itcii-r>v'>y^;^m^Fw«)i^# 

NWcD||Sg^IJ|!cNG{Cj:»JH*tlTV\S, t*6oTII^I 

^gcNG^i^ifli6^^cile^$tlTv^So mtifeflWEGii 

^«^BG^®r^g?fi1I^i:LTStSnTVNS= y 

D-;'?;i/#'UilMGMPEii, xPiJ^^-rsM^^tfT- 

G<Z)tfT'tr*3tlSo 
[0 2 3 6] OTlCfeVNTM^t^gP^I^^^^a-Zt;!/ 

a iB,g,w, P ln,g,wtt^fiKS tlT 1 ocD'^^J' 
c&im^^i^^. mW<zK^V)\yc^t.ir^X(DU-ij 
)l^mmL a 1ml p 1^1 0 1 1 ©g^Ulic 

itm^-t^'^mmiz-^ h u x z * (Z)ttifc-^*tiT v n 

[0 2 3 7] 
5] 

[0 2 3 8] h)\^c&JSiX:^cUtm^ (9 7 a) RXJ^ 
(9 7 b) lCj:'Jtt»$nSo 
[0 2 3 9] cg=G - I (9 7 a) 

cl=T • cg=T • G • I (9 7 b) 

■r3&t)t>vhyi5'xG{Cc}:»j, m¥(D^^^nzMh 
x^^Mi^\H<Difn-jvi^mmm^-^i!b&^t.f)'^^m 

HI iT'XPp{ciU^sti-rv%5yn-/i;i/#ffiil^S: 
^to-Clo®n-5!j;i'#^SJSIR§5:?^fiK't-5>o S^ica- 

fpitstu g^^-^^ltcr^ (9 8) oiii^izim^n 



(22) 
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So 

[0240] ZI=Z* I + Lcl (98) 

fei?)G)ll9^l*ifi^l- ^ (9 3 a, b) ^cJ:S^^- 
7^;l/#^^iJSM7!)^e>^#^4^So #«)Svh'J^'><.3:^^> 

hljU ^ (9 3 b) CD^lC^tl-S^^tC^BS-rSo L 

(9 4 a) fefi^'Wl:^^ (3 4) (D^lZfU^i-^^hLiZ 
J:U#e>tlSo ra)#^lC5t^'J{Cl±^ (3 4) (Z)^5®lt 

[0 2 

[0 2 
[02 



10 



4 1] >l^-^f>X&^-tSXU^>hS:l(f^ 
4 2] ^l^m>U^!7>^T®|■mtS:J^AT^C33V^Tfi 20 



4 3] 
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[0 2 4 4] im^wi.(Dms\HX'iimr':&^. ^^i> 
siiiiE^iB^i:M*y##stis*^e5T'feSo 

StlS. 

[0 2 4 5] fi^TfCfeVNT^StlS^^fiU^liSI 

feVNTttr"= (r", 0", oT) X\ iTJ^^^CJ3^^T 
lir= (r, e, a) ■r*5^$tlTV^So 3llCIHM®t^#. 
Oli'feiS^^' h;br' = (r* , 9' , a' ) ICiU 

^S^I*?IW^T-feSo 
[0 2 4 6] ffim^li. ^^U-:/Big$C (^ (9 0) # 

$>i). z:titcj:'j, ^a-;'?;i/#^ (i=m=o) 

im-t^t:^^(^^ (9 9) *^e>tlSo 
[0 24 7] 
BS5C4 6] 



^ 1=0 m=-l 



(99) 



[0 24 8] 



BS4 7] 



^ (9 9) c±mMo^si)4:^-rt^^^^«±^ m^-iii^^m<nm.mz^y) 



(23) 
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[0 2 4 9] ^oi^^Aet?*6<Jl^«)fe«blc;tf«Jfcidu 

^tb^m'^ii^ (9 9) ©^^^^iCcfcys^aui-* 
[0 2 5 0] mii®^s9#^?:^-rfc«>i-*"^*^^ 

?^^(10 0a), (1 0 0 b) 25.t>* (1 0 0 c) ICJ: 
[0 2 5 1] * 



ii^BIPPll- 3 80 59 

44 

a+=a/ax+i • 3/9y (lOOa) 

a_=a/ax-i • a/3 y (loob) 

az=9/9z (lOOc) 

miC J; ymim^ ( 1 0 1 a) . ( 1 0 1 b) &t>' 

(1 0 1 c) Tb^fiXsi^. 

[0 2 5 2] 



[1SS4 8] 



(lOla) 



(10U>] 



(lOlc) 



[0 2 5 3] tfJit^t>mU^(D^m^ (2M itW^ 85a) , (1 0 2 b) (1 0 2 c) lCj:U^$tvSlO 

mzM^^mft-t^mm^imro^^^Kmrn-^n^^ [0255] 
[02 54] m/^y^^imc^zM Lx mk^ (102 [^4 9 ] 



B.[r^*^ hl^Hfcr)] - •"^"^ cos » 



(102B) 



(103b) 



(102c) 



[0 2 5 6] lll^l^^i:*^^^ (1 0 3 a) , (10 3 [0 2 5 7] 

b) (1,0 3 c) lcJ:y^snsie>lc^®Big5(l- K5C5 0] 



Yf(»,a) 



L+1 



«.m+l 



= a?(21*i)^ 



(S.a) 



(103a) 



(103b) 



(103c) 



[0 2 5 8] r^rdl.±tS^(Dfp(D^k<Dmi (10 4 

O, (1 0 4 b) 2it>' (1 0 4 c) li^fSrefeSo 



[0 2 5 9] 
g5C51] 
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a I fl -t- IB -t- IMI -i -^ 2) 
°x " V (a + 1)(21 + 3) 

1(1 - m -Mjd - m 5^ 

V (a + iMa + 3) 

„ I fl m 4- IMI - m f 1) 
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(104a) 



(104b) 



(104c) 



[0 2 6 0] ±mo>ma^mL<D&m (s s) lo* 5 b) ;s:t>- d o s c) icj:y^$tisii4^K0^^*^ 

#m) *^e5iEM<b^{cS^v^T#e>*iS= ^(102 .^tiSo 
a~c) (1 0 3 a~c) t.^m^-^t>^^:iiL\z [0 2 6 1] 
J: lif?:>:®#®R{C?t>N-S?^ (105a), (10* GgSC52] 

e+[4=^)(Kr) vf (*, «)] =■ (a + 1) «) ^ 

+ kl4^1i(kr) Bin » ei« T^»(», a) ( lOSa) 

a_[4=\icr) vfis, «)] = -t!^(2i * 1) vf;ii(a. «) + 

+ k4^2l(Jcr) Bin » e-l-» Yf(8. a) <l05b) 

e,[4='>M Yf (8, «)] - + a) vJ>Vi(», a) + 

+ khfli(kr)co8«Tff(».a) (105C) 



[0 2 62] ±mDmm^^i^ ( l O e > J: y €. r nrigT'«>So 

tiS^A>^;i/BiKicM-rS!l»^^i5?feM^^i:. [0 2 6 3] 

(10 7 a). (1 0 7 b) (1 0 7 c) tCAU^ [|g5C5 3] 

Bin » ei« YfCB. a) = ajTr^ • ^i^C*. a) - Q? • Yf a) 

(107a) 

Bin » Yf {», a) = ^ • YfTi^*, a) - bf ■ Yi^-i^Cft. a) 

(107b) 

COB « Yf {«. «) " cf.i • Yi».i(». a) - eP • YfVi(». a) (I07c) 



[0 2 6 4] ZniZja^jm^iZltm^ d O S a) , 

(1 0 8 b) 2^0^ (1 0 8 c) iz^^)m^ti^ms'^<m 



[0 2 6 5] 
[^54] 
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(25) 
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[0 2 6 6] 



8_[4^)(lcr) Yf (8, a)] = -J^ * 4+1 'ff+'l (»• «> " 

-l^J'^l«»ifllOcr)vJ'j£^»,a) (108b) 

e,[4^)()cr) Yf (». a)] = -c? Jc 4+1 *r+l(«' «) * 

+ =1-1 •« 4->'^) °' * 

tztbiZ^ (9 9) *^r' .(»•,«' imLXTSmSH-h^kUX*)¥S ( 
1-1) :*sm3 1 :^ff>ffi3!iiaJJ?**«»» ( l + 1 ) :*<io«G(ia^^W«^«U:: 



[0 2 6 7] #mT>'l'^y:^-^ 

[0 2 6 8] BIi 3lcii>r>'T-*^y^'^^^Fw^^l^* 

-A}v^mmmGUPE\t. x^\ii^^-t^mm\^X'nt) 
z.tnzMi.xu-ij)v^mmmitwmmmmBGiH 

[0 2 6 9] ^^im(D±^^it^m.FMMr)VriV:K 

[0 2 7 0] ^^':J-^«yi'>5:#^lciiJ;U:^^vxg|5^I 
[02 7 1] 01 lZltm^^'^^^:^<^^^ *V ^S: 40 



[0tS®fSS#^^^l 



[05] >f >l^-^>:^S:^"rSXl/^>h{CfeW-5, 
[06] m^'WLT.v'V^Z, 3jRtJf4®fc«)(Z)gP4^Mi|)c 

[07] ^-fb-tsxb^^/h-vilStcfcWSSP^^^i^® 

[08] ^J^Sr, e, a $:S^1-|llSi^0T'^S„ 

[09] ^u-A)\^immm(Dtz^(Dmimi'$m.^^ 

[010] m^)\^- 4 c>®S|54>^^(Z)^Mm® 

[011] il^S^S.t>T^gSii^$tiTVNSg|J:J>fUgc?: 
^■tllfSi®0T'feS= 

[012] i^m(Dt:^^(mM^mm.^B^^^x^ 

So 

[013] 2 |gFMM®fefe®iffl^i^^ttJT^^^$ 



(27) 
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[01 1] 



□□□□aanaaa 
□aQQaaaana 

□□□□□ i^!; go 

□□□□□ O: 

aaaaanaaDQ 

□□□□□aDQQD 



□ « WCBG 

H NMM«N6 

[T] ifa^^vu^mmm gmpe 

[T] P-*^I.*«IBB LMPE 



[Ell 3] 



m" if 



□ O : DD 



_ \m 
13 ^n-^\/u9mmm gmpe 

[£1 n-^A'^tUIBLMPE 



